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Abstract 

o 

(N 

The dynamical effects of the QCD-based gluon field overlap factor / (exp{—bskf X]i<, ''?j)) has 

been seen for the scattering process of D° D^ — > u>J/tp. The value used for the parameter kf in 

the exponent of the factor-/ has its origin in numerical lattice simulations. We have calculated the 

cross sections and phase shifts for the scattering process of D°D — )■ uij /ij} when / is used in 

the off-diagonal elements of potential energy, kinetic energy and normalization matrices of a model 

of q^lf systems in the gluonic cluster basis. Our cross sections can be potentially compared with 

i-^ ' experiments. As D°D'^ is considered to be a significant component of the meson X(3872), our work 

*T^. niay have implications for the structure of X(3872) in pointing out the roles of this component 

Oj' of it and other in it; we accordingly comment on meson-meson binding in the D° D„. By using the 

D ' resonating group method, we hope to explore the features of the meson-meson system possibly missed 

r^^ . in Born-approximation based treatments for the almost degenerate states of D° Do and X(3872). 

^ \ 1 Introduction 

<^ . 

OO . In the last few years important progress has been made as far as charmonium spectroscopy is concerned. 

^D ' The many missing states, which were expected in quark model, has been seen. Furthermore some unex- 

^^ \ pected states has also been discovered experimentally which are challenging the quark model seriously. 

t^^ ' These newly born states were given names alphabetically as XYZ [1]. An important member of this 

^— ^ , family is the meson X(3872) which is now generally considered [2]as a mixture of D°D^, D^D_ and 

cc. Any effort to understand it, thus, should understand many quantities composed of combinations of 
its components. A direct lattice study of it would have to calculate many corresponding Wilson loops 
before arriving at any conclusion. A more manageable route could be to make separate models of its 
components, find out their consequences and then combine the models to understand X(3872) or other 
^ ' mesons with a complex structure like that. If separate models use Hamiltonians, then for combining the 

components of AT (3872) the suitable choice perhaps remains a non-relativistic quantum- mechanics-based 
formalism. 

Here we take the first step in this scheme and concentrate on models of D°D^. This is a two quarks 
two antiquarks molecule. We have already worked on meson- meson molecules [3j. But that work did not 
take into account the spin and flavour degrees of freedoms. We now specialize the formalism developed 
to the charm, anticharm quarks along with lighter ones. That is what we have to do study the D°D^ 
system. (Perhaps it is worth mentioning here that there are other meson molecules worth studying as 
well. Rather the strong candidates for meson molecules or tetraquark states belong to light scalars [^ [5] 
[6] as well as some hidden charm resonances e.g. A(3872), r(4260), y(4140) and Z+(4430). Some 
times a common name 'tetraquark' is used for both the types {qqqq and qqqq) because the lattice methods 
do not distinguish between the two types [Q. 

"Scalar mesons are special because they have the same quantum numbers as the quantum numbers 
of vacuum, J^'~^ = 0++ [H]. So in contrast to pseudo-scalar and vector mesons, the scalar mesons have 
non zero vacuum expectation values. They belong to Higgs sector of strong interactions and reflect the 
QCD-vacuum structure." "By taking scalar mesons as qq-systems in the naive constituent quark models, 
it is not possible to reproduce the mass spectrum of the scalar nonet [TOl. (There have been observed 
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15 scalar mesons including /o(400- 1200), /o(1370), /o(1500), /o(1710), fc*(1430), ao(980), ao(1430) [TT] 
instead of 9 below 2 GeV.) Thus whereas SU(3) flavour symmetry works well for reproducing the mass 
spectrum of pseudoscalar (J^*^ = ^) and vector (J^*^ = 1^^) meson nonets in the conventional qq 
picture of mesons, this qq picture is not working for scalar mesons ^12,. The excess of scalar mesons may 
be explained by the introduction of non-conventional mesons with internal structure other than the qq. 
So the possible interpretation [10] of the internal structure may be: l-"They may be four quark states 
with strong diquark correlations" [4] [13] [M] [15] 2- " They may be collective qq excitations" [16] [17] [18] 
[H] [20] [21] 3- "They may be mesonic molecules generated dynamically" [^ [23 ] [24 ] [25] [26] [27 ] 4- 
"They may be glueballs" [28].) 

Coming back to X(3872), there are indications that X(3872), r(4140) and y(4260) are tetra-quark 
or molecular states |8j. Author of ref. [29] identified X(3872) being a tetraquark. Ref. iSOj assume that 
X(3872) is a molecular bound state of neutral charm mesons. Not everyone, though, takes the X to be 
purely molecule or tetraquark. But generally it is accepted [2][3T] that the D°D^ molecule is a significant 
(if not dominant) component of X(3872). (We would like to add that such of mixing of a quark antiquark 
state with a meson- meson molecule is not very uncommon: In ref. [32] it is stated that the " Fock space 
decomposition of the a meson reveals that the dominant component behaves with the increase of A^c, up 
to Nc = 6, as a qqqq then the subdominant gg-like takes over for large Nc up to Nc — 20. The glueball-like 
component stays always at or below the 10% level" . The glue ball states do not mix with qq, qqqq etc. 
in quenched approximation ^33". However full QCD reveals that there is mixing between uu -\- dd stats 
and ss state along with glueball and multiquark channels for any value of J^'~^ of fiavour singlet states. 
So it turns out that for the physical spectrum of scalar mesons there is a mixing of qq, qqqq, glueball etc. 

m) 

For understanding A (3872) we have selected to discuss the D°D^ molecule. This is a pseudoscalar 
vector structure, which is also important in its own because a) PgV (pseudoscalar vector) scattering gives 
excellent information about the spin dependent forces b) PsV scattering phase shifts are experimentally 
accessible [32] . In PgV scattering a D°D system may go to an ujJ/i/j system by making an intermediate 
state X(3872). 

There are many works [27][36][37 [38 31] that treat two-quarks two antiquarks through a sum of 
pair- wise interaction and calculate the meson-meson binding in it. For example, Swanson found [27] that 
a DD system may couple to tuJ/ip system in the one gluon exchange model but, according to him, 
the small interaction is not sufficient for a bound molecular state. Through the present work, we want 
to find implications of non only improving the two-quarks two antiquark modeling beyond the sum of 
two body approach but work out a formalism for meson meson scattering that is not limited to Born 
approximation. In a recent work of us [3], we have mentioned a model that is better than a simple 
sum of two body approach and have worked out its meson-level consequences. We had to use a Born 
approximation there. In the present paper, we use a resonating group formalism to go beyond the Born 
approximation and introduce the spin and fiavour degrees of freedoms. A similar realistic meson-meson 
treatment for lighter quarks has been published earlier [3^. But, in contrast to [39] we now address a 
meson-meson system [D°D^) of current interest. We also give, or plan to give, a much more thorough 
analysis of the meson-meson phase shifts and binding in this system. We include also the cross-section 
that are not in [3^ at all. 

As compared to the sum-of-two body approach, we show in our results section how the inclusion of 
the QCD effect in the form of the gluonic field overlap factor / reduces both the phase shifts and cross 
sections, meaning a reduced interaction. 

We should point out in one respect our present two quarks two antiquarks model is not as close to 
QCD as is the model used in [3]. This is because, for the gluonic field factor / of the model, we use an 
old [40j but computationally more convenient Gaussian form. It is expected to be a good approximation 
because we have checked that the binding energies for a four quark system do not change much if we use 
the Gaussian form of / or more refined surface area form {i\[ of /. 

We think our improvement beyond the Born approximation is an important one. This is because 
we are dealing with states {D°D^), J/ipoo and X(3872) that are close in energy. Energy coincidence is 
the main reason behind a mixing through quark exchange in D°D — > J/ipuu because D°D and 
J/ipuj are effectively mass degenerate. For such almost degenerate states approaches similar to a non- 
degenerate perturbation theory (or Born approximation) fail; see ref. [42) . As our resonating group 
method based treatment is not limited to any perturbation theory (or it has some resemblance with a 
degenerate perturbation theory) we hope to include effects that might have been missed in treatments 



(like Born approximation) that that closer to a non-degenerate perturbation theory. 

In section[2]we have defined the Q^Q^ Hamiltonian and the corresponding wave function with reference 
to spin and flavor basis. There we have also given the form of / that we have used. In sec. [3] we have 
calculated the Potential Energy, Kinetic Energy and Overlap Matrices in the Flavor, Spin and Gluonic 
bases. In sec. |3]we have written the coupled integral equations for Q^Q^ system and solved them for 
transition amplitudes. In sec. [5] we have presented our results and gave conclusions. 

2 The Q^Cf' Hamiltonian and the Wave-function 

In ref. [27, it is given that the Born order scattering amplitude for D''D — > ujj jtl) scattering is 
dominated by the confinement interaction (in the present work we are taking quadratic confinement) 
which is in contrast to light meson scattering (which is dominated by the hyperfine interaction) so 
hyperfine interaction for this process can be neglected which we did here in the present work. 

Using adiabatic approximation we can write the total state vector of a system containing two quarks 
two antiquarks and the gluonic field between them as a sum of product of quarks (Q or Q) position 
dependence function ^r(ri,r2,r3, 1*4) and the gluonic field state \k)g along with flavor and spin bases. 
\k)g is defined as a state which approaches \k)c in the weak coupling limit , with |l)c — |li3l24)c, 
|2)c = |li4l23)c and |3)c = |3i2334)c. In lattice simulations of the corresponding (gluonic) Wilson loops 
it is found that the lowest eigenvalue of the Wilson matrix, that is energy of the lowest state, is always 
the same for both 2x2 and 3x3 matrices provided that |l)g or |2)g has the lowest energy [IT]. The 
later calculations |43) of the tetraquark system were also done with a two level approximation. Taking 
advantage of these observations, we have included in our expansion only two basis states. As in resonating 
group method, ^r(ri,r2, r3,r4) or ^^(R-c, R-fe,yA;, z^) is written as product of known dependence on 
Rc,yfe,Zfe and unknown dependence on Rfe. i.e. ^r(ri,r2,r3,r4) = ^c(Rc)Xfc(R-fc)V'fe(yfe,Zfe)- Here Re 
is the center of mass coordinate of the whole system, Ri is the vector joining the center of mass of the 
clusters (1,3) and (2,4), y^ is the position vector of quark 1 with respect to 3 within the cluster (1,3) 
and Zi is the position vector of quark 2 with respect to 4 within the cluster (2, 4). The same applies to 
R2, y2 and Z2 for the clusters (1, 4) and (2, 3). Similarly we can define R3, y3 and Z3 for the clusters (1, 2) 
and (3, 4). Or we can write them in terms of position vector of the four particles (quarks or antiquarks) 
as follow. 

(ri + rv^ - rv2 - rj) , . 

^' = OTT^ ^^) 

(ri+rr2-rr3-r4) 

^■^ = OTT^ ^^) 

Where r = ^^^, rric is constituent mass of charm quark and m is that of up or down quark. While the 
other definitions of position vectors are 

yi = ri - Tg and Zi == r2 - rj, (3) 

R2 = ^(ri + r^ - ""2 - rj) , yj = ri - r^ and Z2 = r2 - rj (4) 

and 

73 = ri - r2 and Z3 = rg - r^. (5) 

Thus meson meson state vector in the restricted gluonic basis is written as 

2 
|*(ri,r2,r3,r4;5))=^|fc)g|fc)/|fc)3*,(Re)xfe(Rfe)a(yfe)a(zfc). (6) 

fc=i 

The interquark potential in each cluster is taken as simple harmonic oscillator potential (see eq. p2|)). 
So the space part of mesonic wave functions remain Gaussian i.e. 
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Figure 1: Topology 1. 
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C/c(Zfe) 
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(8) 



but with different sizes. k=l,2 correspond to topologies in fig. [T]and fig. [5] respectively, while d^i and 
dk2 are the sizes of mesonic wave functions. For a 3-d simple harmonic oscillator, size square is inversely 
proportional to the square root of reduced mass of quark and antiquark in the meson. So, if 

d — ( Jm ] , d' — { mL^ J and d" — ( mj^e J (we are taking constant of proportionally being 

V ,71 + 771 / V 771+777C / V 177 c + 771 c / 

1, this constant of proportionality eliminates in the following ratios oi d "), then we have 
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(10) 



Here d' = dn = di2 is size of mesonic wave function of a heavy-light quark system corresponding to 
mesons {D°D ) in channel-1, d = ^21 for both light quarks and d" = d22 is correspondingly for both 
heavy quarks for mesons (w and J/tp respectively) in channel-2. 

These Gaussian forms of meson wave functions are, strictly speaking, the wave functions of a quadratic 
confining potential. But, as pointed out in text below Fig. 1 of ref. [44], the overlap of a Gaussian wave 
function and the eigenfunction of the realistic linear plus colombic potential can be made as close as 
99.4% by properly adjusting its parameter d. A realistic value of d mimicking a realistic meson wave 
function depends on the chosen scattering mesons and thus in this present work we actually have chosen 
the realistic value of d for J/^. 

As for the Hamiltonian, for /=1 the total Hamiltonian H of our 4-particle system is taken as ^6] 



H 
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E 



It 

27ni 



^«(r,,)F,.F, 



(11) 



K] 



Where rrii [i — 1,2,3,4) is the mass of ith quark and our kinetic energy operator is non-relativistic. In 
above F has 8 components Fi' = -j-, i' — 1, 2, 3, ..., 8 and F*, = -^, Xi> are Gell-Mann matrices. 



For the pairwise qq potential, we have used a quadratic confinement 



r2. 



= Cri + C with i, j = 1, 2, 3, 4. 



(12) 



As for the within-a-cluster dependence of the wave function, this use of the quadratic potential in place 
of the realistic Coulumbic plus linear may change the full wave function. In the within-a-cluster, this 
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Figure 2: Topology2 



change of wave function is found to result in a change of an overlap integral from 100% to 99.4% only 
provided the parameter d of the wave function is adjusted. It seems that by a proper adjustment the 
parameters of the quadratic (or SHO) model can reasonably simulate a qq or even a q^cf' system. In our 
case, this adjustment of the parameters has been made by a choice of actual scattering mesons in the 
presence of spin and flavour degrees of freedom. 



2.1 Inclusion of Spin and Flavor 

In this chapter we are considering a realistic system D°D by introducing spin and flavour in our 
formalism of treating a four quark system. But for computational convenience we are considering the 
gaussian form of / (see eq. (l24)l ) instead of more accurate minimal area form (see eq. ((23)) ). A q^lf' 
system can have a spin of 0,1 or 2. Here we study D° and D with spin and 1 respectively. So the 
total spin in the incoming channel is 1. But for the outgoing channel we have w and J /i) both with spin 
1. The conservation of spin tells that the spin of the outgoing channel must be 1. Our spin states are 



|1).= I^13^24> 



(13) 



and 



|2). = 11^4%), (14) 

where P represents a pseudo-scalar and V represents a vector meson. For the flavour basis we have 
D° = I cm), D = |cu), ix) = -y=\uu-\- dd) and J /ip = \cc) so the flavour states in channel 1 and 2 are 

\l)s^\cu)\-cu) (15) 

and 

\2)f = ^\uu + dd)\cc) (16) 

respectively. 

Now the spin flavour free resonating group method based integral equations generalize to 

Y, f d'y,d^zkUyk)Ck{zk)f{kUk\g{k I (H-E,) I i),\i)S)fXi{iiiMyiKiizi)^o, (i7) 
1=1 -^ 

for fc = 1, 2. Here H = K + 1^^ + J2i=i "^i- The gluonic states are defined in sec. [2]while spin and flavour 
states are defined in eqs. (|13m6l) above. As operator {H — Ec) is identity in the flavor basis, the overlap 



factors are also the flavor matrix elements of the {H — Ec) operator. Thus the above eq. can be written 



as 



Y, j d\^d\kf{mfUyu)Qk{^k)s{k\g{k\H -E,\i)g\i),xmi)^i{yiKi{^i)^^ (18) 

1=1 •' 

for /c = 1,2. 

Or 

Because of isospin conservation and no annihilation effects each flavour overlap factor is unity. We 
may separate flavour overlap from spin and gluonic overlaps. 

2 

1=1 •' 

According to the (2 dimensional basis) model la of ref. [41], the normalization, potential energy and 
kinetic energy matrices in the corresponding gluonic basis only are 

N^( i\ V V (20) 



3 



/ 1 



V = I . ~ ("13 + "24) 9/(^^12 + -^J- Vi3 - ^24 - «14 - ^23 J | (21) 

^ l/(wi2+W34-«13-«24-«14-«23) "t ("l4 + ^23) ' 



and 



,{k\K\ l)s=Nif)h{J2-^)^(fti- (22) 

This is the modification, through the / factor, to the Hamiltonian. 

2.2 Form of / 

Ref. [41j supports through a comparison with numerical lattice simulations a form / that was earlier |45) 
suggested through a quark-string model extracted from the strong coupling lattice Hamiltonian gauge 
theory. This is 

/ = expi-b,kfS), (23) 

S being the area of minimal surface bounded by external lines joining the position of the two quarks and 
two antiquarks, and bg = O.lSGeV^ is the standard string tension [iSJUT], kf is a dimensionless parameter 
whose value of 0.57 was decided in ref. [JJ by a fit of the simplest two-state area-based model (termed 
model la) to the numerical results for a selection of Q^Q^ geometries. It is shown there |3T] that the 
parameters, including kf, extracted at this SU{2)c lattice simulation with /3 = 2.4 can be used directly 
in, for example, a resonating group calculation of a four quark model as the continuum limit is achieved 
for this value of /3. 

For interpreting the results in terms of the potential for the corresponding single heavy-light meson 
(Qq), a Gaussian form 

f = exp{-b,kfJ2rl) (24) 

of the gluonic filed overlap factor / is used in ref. [JS] for numerical convenience and not the minimal 
area form. But for a particular geometry, the two exponents (the minimal area and the sum of squared 
distances) in these two forms of / are related and thus for a particular geometry a comparison of the 
parameter kf multiplying area and corresponding (different!) kf multiplying sum of squares in eq. (14) 
of ref. [15 is possible. We note that, after correcting for a ratio of 8 between the sum of distance squares 
(including two diagonals) and the area for the square geometry, the colour-number-generated relative 
difference for this geometry is just 5%: the coefficient is 0.075 x 8 = 0.6 multiplying sum of squared 
distances and 0.57 multiplying the minimal area. But, as the precise form of / is still under development 
(the latest work [43] has covered only a very limited selection of the positions of tetraquark constituents) 
and the expression for the area in its exponent needs improvement, it is not sure precisely what value of 



the kf best simulates QCD. Furthermore we have checked that the binding energies of four-quark system 
do not vary much at least for the square geometry. So for the present work we are taking / reported in 
eq. (14) of ref. [48] and is numerically easier to deal with. This sec. I2.2l has been taken from our previous 
work [3]. 

3 Potential Energy, Kinetic Energy and Overlap Matrices in the 
flavor, spin and gluonic basis 

We define in this section the potential energy, kinetic energy and the overlap matrices of potential 
energy, kinetic energy and identity operators in the flavor spin and gluonic basis with out doing any 
space differentiations and integrations. First we discuss the matrix elements of potential operator 

4 
V^iqmqsq^) = Y. F-Fj^' (25) 

i<j 

, where Vij is defined in eq. (I12p . in the gluonic-spin bases. 

Ks,i,s = gsik I V'' I l)gs = ^F,.F,(i;,j)fe,.,;,, = ^ a., (/c, /)(«., )fe,,,,Afc, (26) 

Here 0^3(1,1) = a24(l' 1) = ^|' 012(1,!) = 023(1^ 1) = '^2li^^'^) = «34(1' 1) = ^^^d 012(1,2) = 
034(1,2) = —0-^3(1,2) — —024(1, 2), —023(1, 2) = — 0]^4(1,2) — |/ and similarly others are the matrix 
elements of the Fj.Fj operator in the gluonic basis {|l)g, |2)g}. Aki are the elements of spin overlap 
matrix A. 

(See appendix-a for details). 

The hyperfine interactions are neglected here because they do not affect the scattering amplitude 
much [2^. So the potential energy matrix becomes 

yP ^ ( ^K3+^24)ls,ls |/(«12+^^_-f'l3-f'24-f'l4-f'23)l^.2s ^ (28) 

V t/(wi2 + %4 - Wl3 - «24 - ""14 - ""23)2^, Is -Y (VlJ + W23)2s,2s J 

Where the subscript s indicates the spin states. 

Now we come towards kinetic energy operator. In spin space, K and (X]i=i fTT-i — Ec) are unit operators. 
So their matrix elements in spin-gluonic bases are equal to matrix elements in gluonic basis multiplied 
by spin overlaps, i.e. 

Kk,.,i,^ = gs{k I K I /)g, = s{k I Os g{k\K\ I), = Au ,{k\K\ I), (29) 



and 



,(/c I (^m, - E,) I l)gs = {Y^m, - E,) g,{k \ l)^^ = (^ m, - E,)Nk^^,i^^ 

i—l i—1 i—1 

4 4 

= CYm,-Ec) s{k\l)s g{k\l)g^CYm^-Ec)Akl g{k\l)g (30) 



i=l 



where 



Nk,a.s=ss{k\l)gs (31) 

Using eqs. (|3T]), ^ and ^ in eq. dH]), we get 

/'dX[K.,„fe,.(Rfe,R'fc) + Vr,^,,fc^JRfe,R;,) + (^m, -i?,)Nfc,^,fc^^(Rfc,R',)]xfc(RL) + 

i—1 

f d^Ri \Kk^^,^^ (Rfc, Ri) + Vf , (Rfc, R;) + ( V m, - Se)Nfe^,,,^, (R^, R/)! Xii^i) = 0, (32) 



for kj = 1,2. 

Kfeg,,i,,, ^k,,,i^, and Nk^,,!^, are defined as 

I d'R[Kk^^,iJTlk,-El'i)xii'R-'i) ^ I d'ykd'zkUyk)Ckizk)Kk^^,i,M'RiMyi)Cii^i) (33) 

J d^B:iV^^^^iJRk.R[)xi{R[) = J d'y,d'zkUy,)a{zk)V,l^i^M^iMyi)Ci{zO (34) 

y d3R;Nfc^^,i^,(Rfe,R;)xz(R0 = y dVfcd3zfcefc(yfc)a(zfc)iVfc,„i,.»(R06(yi)0(zi) (35) 

3.1 Space Differentiations and Integrations 

In this section we describe how we performed spatial integrations on right hand side of eqs. p3ll35p . For 
that we have to substitute Kk^^j^^, V^ ; and Nk^^jg^ expressions from eqs. (f^ . (^5)) and ([51]) . As for 
as Kk ^j ^ is concerned we have to do spacial differentiations as well because it involves kinetic energy 
operator. In the case of fc = Z, i?; is linearly independent of y^. and z^ thus xiC^i) can be taken out side 
of integration. For k ^ Iwe have to write Rj, y; and z; in terms of R^, y^, and z^. 

First for k = I in eq. (|5T]) by using eqs. ([27] and ^U^ for spin and colour overlap matrix diagonal 
elements, we have 

^fcg.,fcg, = gs{k I k)gs = s{k I k)s g{k I k)g = ^fefe N{J)kk = 1- 
Putting this in eq. (|55|) we get 

J d^R',N,g^,,gSRk,K)XkiK) = XkiRk) J d^yJ^ZkaiYkKlM. (36) 

Because £,kiyk) and Cfc(zfc) are normalized, so we get 

I d'R',Nkg^,k,^(Rk,R',)xki-R.'k) - Xfc(Rfc). (37) 

Or 

Nkg^,kJRk,K) = mk,K), (38) 

for fc = l,2. 

Now for kinetic energy matrix elements from eq. (1291) for k = I, we have 

Kug,,k,,^Auk g{k\K\k)g. (39) 

■Akk ~ 1 from eq. (P7)) while from eq. (P^ we have 

4 9 

g{k\K\ k}g^N{f)l,[j2-§r)^(fik (40) 

1=1 ' 

Using eq. ([20]) 

.(^i^i^). = (E-& (41) 



As the four quarks are not of the same mass in this chapter so there appears index i with m. 
Substituting 






^ v^ 






i^-m '^^) 



for k=l,2, in eq. ([33)) we get 



4 2 

dVferf'zfcefc(yfc)a(zfc)(^-^)xfe(Rfe)6(yfc)a(zfc) 

i—1 * 



(43) 



By using eqs. ([T|), ([5]), ([3]) to ([5]) , for Ri , R2 , R3 and for y^, and z^, we can write in the set {R^:, y^., z^} 

1 



(E 



* V? 



2m, 



2m 



kY^,+<ZfcY^ +tfcV^J 



(44) 



for fc = 1,2, 

with 771 being the constituent mass of the hght quark, up or down. 

Here 

2 r+1 r+l 2 

si = — — r, <7i=ti = , S2 = — — , (72=2, t2 = -. 

r + l r 2r r 

By using eq. (|44p in to eq. (l43l) and doing all the space differentiations and integrations, we get 



Kfc,^,,fc^^(Rfe,Rl,) = J(Rfc,R',)g(A,i + Afe2) - ^Y^, 



for fc = 1,2. Here 



Afci — 



gfe 



2md2^ 



and A 



tk 



fe2 



2mdl^ ' 



For the potential energy matrix, from eq. ([26| and eq. (|27|) for fc = Z, we have 



(45) 



(46) 



(47) 



(48) 



i<j 



The only non zero coefficients aij{k,k) for fc = Z = 1 are 013 and aj^j and for k = I = 2 are 014 and 023. 
Using the quadratic form of potential from eq. ([12]) , eq. (|48)) becomes 



Using this in eq. ([34l) we get the result 



2C + Cy2 + Czt 



vi 



S{Rk,K)[-^C-4C[dl,+dl,] 



(49) 



(50) 



For writing the off-diagonal elements firs we write y; , z/ and R/ as linear combination of y^, , Zfe and 



Rfe. This means, we use 



R 



■2 — 



Yi-zi 



R3 = -(^1R 



Ri = 



Ri 



2 


2r 

1 


' (r + iy 


ya - rz2 


r + l 


y2 + ^Z2 


r + l 


2r 

^ + T ^-r 



(yi + zi) 



R3 



;^)R3+(^f^(y3 + ^3) 



(51) 
(52) 
(53) 
(54) 
(55) 



R2 = y^ (56) 

It is clear from the above relations that the integration variables y^. and Zfc are not appropriate for the 
integrations to do next. So, we have to replace them by their combinations with one being identical to 
R; and other independent of it. Because of different quark masses now we can not write e.g. R3 simply 
as linear combination of y^ and Zi. Actually in place of Ri, y^ and Zi, we used 

fr^R3^-(^)R.+ (^(y.+zO (57) 

r+lJ ^ (r + iy 



gi = -(:^^)i^i-7^^^(yi + ^i) (58) 



R2, y2 and Z2 are replaced by 



hi = R2 = ^^^-^ . (59) 



f2 = R2 (60) 



y2 + rz2 
§2 = ^3= ^^^ (61) 



h2=Ri = ^^^-^. (62) 

r + 1 



Similarly for R3, y^ and Z3 we have 



f3 = Ri = -(^)R3+(;:|^(y3 + z3) (63) 

g3 = -(:^)R3-(^(y3 + ^3) (64) 

ha = R2 = ^^2^ . (65) 

With these new integration variables f/c , g^, and h^ (or their combinations ffc + g;, and f a; — g;, ) for k — 1,2 
we can write expression for N^ ^./ ^ and other kernels present on L.H.S. of eqs. p3ll35p . First of all we 
write expression for Nk ^.i ^ from eq. (|3ip with the help of eqs. ((27l and [20]) 



^fe,.,i,. = ^Aki (66) 

for k ^ I. Substituting in eq. p5p expressing all the vectors there in terms of ffc, g^, and h^ and doing 
all the integrations other than R/ we get 

Ni,.,2,. = ^exp{-hRl - hRl) (67) 

and 

N2 1 =Ni 2 (68) 

Here 

lo = ir + l)ir^2i{Tra(f)^ (69) 
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and 



1 /r + l\2 



a 



h=Ak 



1 / 2r 



2d2 V r + 1 



a = 8kd^ 



r' + l 



+ l + r — 



(r + 1)^ 



(70) 
(71) 

(72) 



Sfed-" 



r^ -1 



l + r" 



r^-l 



(73) 



7 = Sfcd^ 



H + 1 



+ l + r" 



(r-l)2 



(74) 



.V2F+1) 

The above integrations leave the kernels as functions of ft, g;^. and hfc which are then replaced by Rfc 
and Ri using eqs. (|57ll65p . E.g. for Nig^^2g^, fi — gi is integrated out leaving with 

fi+gi=-2(^')Ri 

and 



etc. 



.r + 1 
hi=R2 



Now for the kinetic energy kernel for k ^ I. Here along with integration firstly we have to do space 
differentiations as well. From eq. (P^ using eq. (I?fl) and eq. (|22p . we have 



-.^....^j^/HE-fi)/'- 



(75) 



We have expressed ( X]i=i ~2m" ) ^^ ^^- ©1) in term of V^^^, Vy and V^^. Now using relations from 
eqs. (|5m65p . we can write for fc = 1 

^ — ' 0-r-n I 0-r-n 



1=1 



8(r-^ 2 

U _(_ 1)3 — (fi+Si) 



32r ^2 



r + 1. 



(r + l)3-(^i-si)^ 2r 



-vL 



(76) 



and for k — 1 



V-'-^ 2mi/ 2m 



r + 1, 
^7 



"-^f2 



(r + l)2 



Il-^(g2+h2) T'JL(g2-h 



rV; 



Or 



(E 



' V? 



1 



Irrii) 2m 






r + 1 

~^^R.2 + (r + 1)2 i^(R3+Ri) + ''^(R3-Ri: 



Substituting eqs. (|77)) and ([75]) in eq. ([33]) . we get 

L 1 



and 



Ki 



K2 1 



2m3V2 



/o 1 



riiRf + ri2R2 + no exp{-hRi - ?2R2) 



2m3V2 



r2iRi +r22R2 + r2a 



expi~hRl ~ hUl), 



(77) 



(78) 



(79) 



11 



where 



rii 



1\4 8(r- 1)2 f /r-l\/ 8fc I + ^/^ \ { P r-l> 



V 2 / (r + l)3 |Vr + l/V(r-l)2 ^ (r- 1)2^2 J " Va " r + J 

/2fc 1 \1 32r f/r-l\/2fc 1 \ 

W "^ 2dVRTTT7^/ I ^ (r + l)3l Vr+JVT"'' 2d2,/r(l + ^)j 



r 2dPy/r{l + y/^)J \ (r + 1 

\a r+l/V 



2 



/3 r-lNj^-r2 + l ^-3/2^;^ 
~^ + 4d2 



^'lO 



3 / r + 1 
2 V 2 



,(--1)^ 



1 + v^ 



32r 



(r + l)3 V(r-1)2 (r- 1)2^2; (r + l)^ 
2 f ai'^ i\2 /oi; I ^2 



r2 4^2 I [ ^2 a^ M (r + l)3 V ?• 2^27^(1 + 7^) 



32r / fc(r2 + 1) r" + 1 



(r + l)3 



4d2 



6(r + l) 
2^^ 



2A: 



1 / 2r 



2^2 V r + 1 



(80) 



(81) 



and 



ri2 =4 



r + 1 
2r 



2fc 



1 / 2r 



2^2 V r + 1 



?'22 = ?'12, 



1\2 /■ k r 2 \2 



r20 = 



(r + 1) 
r + l\2 



1)2V 2 / a 1 V 4^2^ Vr2 4^2^ f "^.--^^2 



- 1 \ / fc r 2 

''+4d^J+''V^+ 4d^ 



6^ 2fc 
2r I 



(r + l)2 



1 / 2r 



2d2 V r + 1 



1 \2 



/?\2/fc r2\2 



--<-')1('-!^)(^-i^)-(-9(^-^ 



(82) 
(83) 



(84) 



(85) 



Now for the potential energy kernel usmg eq. (pS)) . eq. ([T^ . and eq. ([27]) in eq. (12^ . we get for A; 7^ Z 

(86) 






2,2 2 2 2 2 

Ya + Z3 - Yi - zi - y2 - Z2 



and 



< 1 = < 2 



(87) 



Putting this expression in eq. p4p . changing variables and doing all the integrations, we get 



Vr , = — CN 

i^gs -.t-gs 



Q ' i^g s ;'' gs "''' f^gs t'' g s 



where 

Vig^^2g, = V2g,ag, = lo niKl + no expi-liHJ - hltll). 
^kg,,ig, is defined in eq. (1571) and eq. (p5|) . Here 

8 /r + lx2d2 
6 \ r / a 
and 



(89) 



(90) 



ni 



4^f (r + l)4 W/3 r-1 
9 I 4r2 /la r + 1 



(91) 
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4 Writing Coupled Equations 

For fc = 1 putting expressions from eqs. (gni, 1^, 1^, ((75)) . 



and ([57)1 in eq. ([5^ . we get 



^^(Aii + A12) - |i^V^^ ~^C- iC[dl, + dl^] + 2m(r + 1) - Sj] xi(Ri) 



^2^0 / d^R; 



1 1 

2™ 6 



riiRj + ri2R2 + rio + - niRf + Uo 



-i(^e + ^C-2m(r + l)) 



ea;p(-;iR?-;2R2)X2(R2) = 



(92) 



Now for fc = 2 putting expressions from eqs. gS]), dSO]), dMl), dlSl), dMl), jMl) and §8^ in ([32]), we get 



^7(A2i + A22) - ^^l, -^C- 4C[dl, + dl^] + 2m(r + 1) - E,] X2(R2) 
4 2m 3 



y2Zo / d^Ri 



1 1 
2to6 



r2lRi + r22R2 + f20 



1 r 
2 



niRi + rio 



-K 



8 



E, + ^C-2m{r + l)^ 



ea;pHiR2-/2R2)xi(Ri) = 



(93) 



It is clear from eqs. (|92l) and ([M)) that off-diagonal parts vanishes for large values of Ri and R2. So 
in this limit the two equations become 



'^Aii-^C-4Cd?i+m(r + i; 



+ (^Ai2 -\c- ACdl^ + m{r + 1] 



-(|j;^.+^' 



Xl(Rl) ^_Ri^oo 



(94) 



and 



V4 



A21 - ^C - 4Cdli + 2m 



(^Xll 



^k, + E. 



+ (?A22-xC-4Crf2 



.4-'^^ 3- 



22 



2m, 



\2m' 

X2(R2) >R2 



(95) 



For the large separation there is no interaction between the two mesons. So the total center of mass 
energy in the large separation limit will be the sum of kinetic energy of the relative motion of mesons 
and masses of the two mesons, i.e. for i?i — > 00, we have 



1 



2^ 



and for R2 ^ 00, we have 



D°D 



1 



V^, + Md + A%o 



-Y.R, + M^ + Mj/^ 



xi(Ri) = £;cXi(Ri), 



X2(R2) = ^cX2(R2)- 



By making a comparison of cq. ([M]) and eq. ([M|) . we can write eq. ([M)) as follow 
1 



(96) 



(97) 



Md + A%o. 



^f^D°D'' 



1 r 
'2 



- V^, - Ec 



1 



X 



l(Rl)+V2/o / d^R; 



1 1 

2m 6 



riiR^ +ri2R2 + no 



niR^ + Uo - -{^Ec + -C - 2m{r + 1) 



expi-hnl - ;2R2)X2(R2) = 0. 



(98) 
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Similarly by making a comparison of eq. (1951) and eq. ((97)) . we can write eq. ((93)) as follows 
1 



M^ + Mj,^ - 



J/^ 



S^wJ/i/- 



V^. - Ec 



X2 



(Ra) + V2I0 / d=^Ri 



1 1 
2m 6 



r2lRi +r22R2 +7-20 



+ : 



1 r 



niR^ + Uo 



\{Ec + \c-2m{r + l) 



expi^WHl - /2R2)xi(Ri) = 0. 



By taking Fourier transform of eq. (IMl) and (IM)) . we get 

Mz5 + i\%o. 



and 



Here 



1 ^12 



^ -p?-^;. 



2m 6 



■B2{l2 



Xi(Pi)- 

J_!i£ + !!^ _ :^ 
2m 6 2 6 



2m 6 2 



A2a2)^b(Pi,/i) 

^2(^2) F„(Pi,/i)=0, 



M^ + M, 



J/4' 



2MwJ/v 



P2 -Ec 



X2{P2) - T^'^A,{h)Ft{P2j2) 

Zm D 



2m 6 2 6 



Aiih 



2m 6 2 



r)siGi)li^a(P2,/2)=o. 



(i'^Ri 

Xfe(Pfc) = / ,„ ,3 exp[tPfc.Rfc]xfc(Rfc)> 
(27r)2 



(99) 



(100) 



(101) 



(102) 



Where 



E'^^E, + -C-2m{r + l), 



Ak{u) = C / d-^Rkexp[-uIli]xki'Rk), 



I' = 72^0, 



Bk{u) = C I d^Rkexp[^uRl]Rlxk{Rk), 



(103) 
(104) 



(105) 



(i'^Rt 
Fa(Pfc,u)= / ——-jexp[LPk-Rk]exp[-uRl], 
[2tt)2 



(106) 



d"^Rt 
Fb(Pfe,M)= / -—-^exp[t.Pk.Rk]Rlexp[~uRl], 
(27r) 2 



(107) 



and 



for fc = 1,2. 



FaiPk,u) = 3 exp 

(27r)2 



4u 



Fb{Pk,u) = Fa{Pk,u) 



1 

2^ 



3- 



2u 



(108) 



(109) 
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For the incoming waves in the first channel eq. (|100p and (|10ip can be formally solved [40] as 



and 



where 



XiiPi) 



p'cW 



Ai(pi) 



2m 6 2 



FaiPuh) 



Ai(pi 



2m 6 ^^'^'^ \ 2m 6 2 6 ^ ^ '^ 



^2(^2) 
= 



, . Fh{p2,l2) 1 r22 . ,, . Fa{p2,l2)\( 1 7'20 , «„ E^^ 



A2(p2) 2m 6 



A2(P2 



2m 6 2 6 



AiGi) 



2m 6 2 y 



)Bi{h) 



Aibi) = 



p? 



^l^D°D° 



Md + Mj^o, -Ec-LC , 



(110) 



(111) 



(112) 



A2(P2) = 



pI 



2^*^,7/1/. 



M^ + Mj/^ -Ec-ie 



(113) 



and 



p,{l) - ^/2/i^„^o.(S, - Mz3 - il%o.) 



Pc{2) = ^2^i^j/4E,^M^-Mj/^). 



(114) 



(115) 



Because of the spherical symmetry of the S-wave (/ — 0) here, P^ is replaced with pi (magnitude) with 
i = 1, 2. Using the Parseval relation corresponding to eqs. (|104p and (jlOSp . we have 



Akiu) = AttI'^ / dpkplFa{pk,u)xkiPk) 



and 



Bfc(u) =47r?^ / dpkplFbipk,u)xk{Pk), 



(116) 



(117) 



for fc = 1,2. 

Multiplying eq. (jllOp by 4:TrplFa{pi,li) and integrating w.r.t. pi and with using eq. (|116p we get 



Alii 



^ = 4^F„(p,(l),/i) - A2{l2)wil^ - B2{l2)W^ 



(1) 

2 



(118) 



where 



M^r^J^ 



2m 6 2 



J((^i,/i) 



J_!!l2. _|_ !!£ _ :^ 
2m 6 2 6 



Ji(Zi, Zi 



(119) 



W, 



(1) 

12 



1 ri2 
2m 6 



'^i(^ii 'i)j 



(120) 



and 



Jk{u,v) =Att I dpkPk 



2Fa{pk,u)Fa{pk,v) 



^k{Pk) 



4(^*7 «) =47r / dpfePfc 



2Fa{pk,u)Fh{pk,v) 



^k{pk 



(121) 



(122) 
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Similarly multiplying eq. (IllOp by 47rp^i^b(pi, ?i) and integrating w.r.t. pi and using eq. (|117p . we get 

^^ = AnFt{pc{l),h) - A2{h)W^l^ - B2{l2)W^l\ (123) 

where 



<) ^ 



2to 6 2 



j{'ai,;i) + 



2m 6 2 6 



J{(Zi,Zi) 



(124) 



and 



W. 



(1) _ 1 ^12 w/; , ^ 



J'ki'u-iv) =47r / dpfePfc 



2Fb{pk,u)Fb{pk,v) 



Afc(pfe) 



(125) 



(126) 



In the same way multiplying eq. (jllip by ATip^Fa{p2: h) and 47rp2^(,(p2, ^2) and integrating w.r.t. p2 and 
using eqs. (|116l and lll7p . we get 



^2(^2 



.(2) 



Ai(?i)W^i\^^-Bi(/i)W^i^2^ 



(2) 



(127) 



and 



B2{l2) 



.(2) 



^l('l)W^21 -5l(^l)W^22 , 



(2) 



(128) 



where 



•<-i1^'''<'-w + 



2m 6 2 6 



^2(^2,^2 



(129) 



and 



Wi 



(2) 



2m 6 2 



J2(/2,/2 



W^il''=-7^?^2'('2,/2) + 



2?77, 6 



2m 6 2 6 



^2(^2,^2) 



<^) = 



1 r2i 
2m, 6 



ni 
2 



>^2('2,'2)- 

Equations piSp . p23p . p27p and ()128p can be written in the matrix form as follows 

WVi = V2, 



with 



W 



W[^ \ 



w. 



(1) 



/ C 

v;' 

V w^f wf^ v:' j 






22 





(130) 
(131) 

(132) 

(133) 



(134) 



Vi 



( A,{h) \ 
Biih) 

A2il2) 
\ Mh) J 



(135) 
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and 



Vo =47r 



/ Fa{pc(l)M) \ 

Fb{pc{l)M) 



V ; 



From eq. ()133p . we can have 



Vi = W-^V2, 



which gives values of v4i(^i), i?i(^i), ^2(^2) and -62(^2)- So cq. (|110l) and (|llip become 



and 



where 



XI (Pi) = 



5{p^~p,{l)) 1 



P?(l) Ai(pi) 



T^l^'^A2(/2) - W!^^^B2{12) 



w['^ = 



X2{P2) 



2m 6 2 



A2(P2) 



A2) 



W^'>Ai{h)-W^'>B,ih 



(2), 



Fb(pc(l),;i) + 



2m 6 2 6 



(136) 



(137) 



(138) 



(139) 



Fa{pc{l),ll), (140) 



i^i^^ = -^^F.(p.(i),M, 



(141) 



and 



(2)^_ 1 ^^^^^(2),?2) 

2m 6 



2m 6 2 6 



J^a(Pc(2),/2) 



w. 



(2) _ 



2m 6 2 



J^a(Pc(2),?2). 



(142) 



(143) 



From eqs. (|138p and (I139P we can read off the T-matrix elements Tn and T21 [10] ■ These are proportional 
to the co-efScienl 
and T-matrix is 

s = 1 - 2tr 

when both the channels are open. Eq. (|144p gives 



to the co-efficient of Green's function operators — -^-t — 7 and — -^-t — y Here the relation between s-matrix 



(144) 



and 



where, 



and 



Tn = 2fij,„j^o,^p,{l)[wi'^A2{l2) + W^^^B2{l2] 



T21 = 2^^j^ Jpe(l) J- Kf ^Ai(/i) + W^i''Bi(/i) 

^ V ^1 '- 



Wl 



1^2 



Pcil) 


^^D°D'" 


Pc{2) 



^^LoJ^p 



(145) 



(146) 



Similarly T22 and T12 can be found, using 
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V-2 =47r 



( ' \ 



Fa{Pc{'2)M) 

\ Fb{pc{2),h) J 



(147) 



in eq. (|133p . for incoming waves in tlic 2nd channel. These become 



.(2) 



.(2), 



T22 = 2^l^J^-p,{2) W^rMi(/i) + W^^'B.ih] 



and 



TT 



T12 = 2^^„^o.^p,(2)./^ W^i^'^A2(l2) + M/i'^i?2a2) 



(148) 



(149) 



4.1 Fitting the parameters Used 

By making a comparison of eq. ()97p with eq. (|95p . we have 



M 



J/i' 



rAo 



C - 4Cd^2 + 2mc 



(150) 



From — 4C/3 = /XccA52/2 . Also reduced mass for a pair of charm and anti-charm quarks /Ucc — tiIc/^ 
because c-quark and c-quark have the same constitutional mass. Thus — 4C/3 = TOcA|2/4, along with 
A22 = l/mcd22 gives 



4C/3 = A22/4d^2 °r - 4C<i^2 ^ ^^22 1 ^ 



(151) 



We have fixed the parameters uj and C by spin averaging of masses of charmonium in the state IS" 
and the state 25. 

From eqs. (|150l and ll5ip . we have 



3Afj/^(15) + Af,,(15) _3^ ^^, 
^ 2^22 + C 



and 



(152) 



3M^(2g) + M^^(25) _ 7 



A22 + C 



(153) 



The factor 7/2 with A22 in above equation comes from the relation 
for n = 1 and I = 0. Whereas in eq. (|152p above, the factor 3/2 with A22 



Where, C" = -|C + 2™, 
i^„/m = A22(4n + 2Z + 3)/2[49 
is for ri = and I = 0. 

Put the values of masses Mj/^(15') = 3.09691606^, M^^(15') = 2.9803Gey, M^{2S) = 3.6860906^ 
and M^^(2S') = 3. 63706 V" from (PDG) ref. |9i in eqs. (|152p and (|153p and solving them simultaneously, 
we get A22 = 0.3G3G28GeF and C" = 2.61322Gey. We took constituent mass of light quark as m = 
0.33GeF and that of charm quark m^ = lA7MGeV [50]. This gives G = 0.259185GeF. The relation 
A22 = l/mcd22 then gives the value of ^22 or d" = 0.294583fm and using this value in eq. (fTOj) gives 
the value of d2i = d = 0.428648fm and by putting this value of d in eq. ([9]) we get du = di2 — d' — 
0.379058fm. By putting the above values of ^22 and A22 in eq. (|15ip we get the value of the parameter 
G = -0.0254714Gey3. The value of 5, is taken as O.imeV^ [47]. We take the value of kf = 0.075 [48]. 
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4.2 Finding the Phase Shifts, Cross Sections and Binding 

To find out the phase shifts we have used the same relation 

s = I~2iT (154) 

and [51] 

s = exp(2iA) (155) 

and by comparing their right hand sides, we have 

1 \ ^.f Til Ti2\ f 1 \ , ^. f Sn 5i2 



i)-"{n\ r;:j-(o \}^^'{z z) "^«' 

6ij — —Tij for i,i = 1,2 |51) . Because the phase shifts are much less than a radian, so we can neglect 
the higher powers of exponential series for eq. (|155l) above. To find out the total cross sections we have 
used eqs. [27 to 29] of ref. [52] while according to our definition of Fourier transform hji has been taken 
equal to Tji [53] [51] for z, / = 1,2. The total cross-section is 47r times differential cross section for any 
value of centre of mass kinetic energy, T^. The condition for the existence of a bound state of the q^<f 
system is [33] 

detW^ = 0. (157) 

For the definition of matrix W , see eq. (|134p . 

As this is a determinant of a 4 x 4 matrix, once we select the numerical values of the parameters in 
it, it would have more than one values of the energy that satisfies it. 

5 Results and Conclusion 

Fig. [3] represents the elastic scattering phase shifts in degrees against centre of mass kinetic energies in 
GeV for the process D°D — > D°D , corresponding to t- matrix element given in eq. (|145[) . Fig. [5] 
represents the elastic scattering phase shifts for the process ojj /^ — > uJ /tJ), corresponding to t-matrix 
element given in eq. (|148|) . In fig. [4] we have shown the inelastic scattering phase shifts for the process 
D°D — > ujj lt\), corresponding to t-matrix element given in eq. (I146p . The same results are obtained 
for the inelastic scattering process ujJ /ip — > D°D , corresponding to t-matrix element given in eq. 
p49p . The reason for selecting w meson in the scattering process is that the mass sum of w and J /ip is 
almost equal to the mass of X(3872) meson. To find out these phase shifts we have used the relations, 
between s-matrix elements and corresponding phase shifts, given in the reference [5T] . 

We have found a bound state (see eq. (|157p . condition for the existence of a bound state) at total 
centre of mass energy Ec = 3.6942Gey. But we can get a bound state at Ec — 2>.S72GeV by multiplying 
the off-diagonal terms of potential energy, kinetic energy and normalization matrices (i.e. multiplying lo 
of eqs. (j^ andlM )) ) by a factor- Yq for value 1.97. See fig. [S] where we have ploted the factor Yq against 
total centre of mass energy for a possible bound state. (As the equation, last in the above section, that 
gives this energy value can have more than ne solutions there may be other curves for factors like this 
Yo as well.) This bound state at E^ = 3.872Gey may be recognized as X(3872). In fig. |3]we have -ve 
phase shifts indicating an attractive potential so the possibility of D°D bound state being X(3872) is 
very much alive. Alternatively if we take the value of fc/ being 0.0048 instead of 0.075 then we can get 
a bound state even at 3.87206^. So by adjusting the free parameter fc/ we can get the bound state at 
2,m2GeV. 

In figs. ((71 [8] and [9]) we have shown the total cross sections for the scattering of D°D — > D°D , 
D°D — > ujj /-if), and (joJ/ij] — > w J/i/', if the same eqs. [27 to 29] of ref. [55] are used for the inelastic 
process as well. These results can be compared with the experiments. We have also calculated the total 
cross sections for D°D — > J /ijjtt" and vise versa scattering. 

The same technique may be used to study the other scalar mesons e.g. Ar(3940), F(4260) and Z(4433). 
In future we plan to study the structure of scalar mesons by using the more realistic linear-plus-columbic 
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Tc(GeV) 

Figure 3: Elastic scattering phase shifts for the process D°D^ — > D°D^ 



potential and corresponding wave functions along with the use of the surface area model of / of eqi23l 
and not Gaussian / of eq. [24l In figs. (fTTJto lTSl) phase shifts are shown for kf = 0, while in figs. (|14lto 
[TB)) total cross section are given for kf — 0. It is to be noted that both phase shifts and the cross sections 
decrease when we use the QCD effects in the form of /. 
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Figure 4: Inelastic scattering phase shifts for the process D°D — > ujj/^p. We checked that the same 
results are obtained for the inelastic scattering process utJ/tp — > D°D . 
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Figure 5: Elastic scattering phase shifts for the process ujj/tp — > coJ/ip 
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Figure 7: Total cross sections for the scattering process D°D^ — > D°D^ 
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Figure 8: Total cross sections for the scattering process D°D — > coJ/ip, if the same eqs. [27 to 29] of 
ref. [52] are used for the inelastic process as well. We checked that the same results are obtained for the 
inelastic scattering process ujJ/^p — > D°D . 
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Figure 10: Total cross sections for the scattering process D"D — > J/iJ^t^, if the same eqs. [27 to 29] 
of ref. [52] are used for the inelastic process as well. The same results are obtained for the scattering 
process J/tpn — > D°D . 
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Figure 11: Elastic scattering phase shifts for the process D°D^ — > D°D^, for kf — 0. 
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Figure 12: Inelastic scattering phase shifts for the process D°D — > cuJ/'ip. We checked that the same 
rcsuhs arc obtained for the inelastic scattering process ujJ/ip — > D°D , for kf = 0. 
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Figure 13: Elastic scattering phase shifts for the process ujj/tp — > ujj/ip, for kf ~ 0. 
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Figure 14: Total cross sections for the scattering process D°D^ — > D°D^, for fc/ — 0. 
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Figure 15: Total cross sections for the scattering process D°D — > ^J/'ip, if the same eqs. [27 to 29] of 
ref. [52] are used for the inelastic process as well. We checked that the same results are obtained for the 
inelastic scattering process utJ/tp — > D°D , for kf ~ 0. 
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Appendix A: The spin basis 

For total spin 1 different possible spin states for a meson meson molecule are 



I%^24>-^(|1,1) + |1:0) + |1,-1)), 

V stands for a vector meson. Inserting the completeness expression for an identity / the above equation 
may be written as 

1^13^24) = ^(|1, 1,1,0)(1, 1,1,01 + [1,0, 1,1)(1,0,1,1|)|1,1) + -^([1,-1,1,1) 

(1,-1,1,1| + |1,1,1,-1)(1,1,1,-1| + |1,0,1,0)(1,0,1,0|)|1,0) 
+ -i=(|l,0,l,-l)(l,0,l,-l| + |l,-l,l,0)(l,-l,l,0|)|l,-l). 



Using Clebsch-Gordan coefficients it becomes 

|1,0)|1,-1)-|1,-1)|1,0)) 



l%V^24>-— (|1,1>|1,0)-|1,0)|1,1) + |1,1)|1,-1)-|1,-1)|1,1) + 
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Again using identity 



/ = 



and using the Clebsch-Gordan coefficients we get 

i%F24) - -^ [ tm + nn - mt - im +V2 nn -V2 ;nt + m; 

V 1^ '- 

+ un - ;n; - ;;;t 

From eq. (I158p by interchanging positions of 3 and 4, we have 



1 1 1 l\ /l 1 1 1 


+ 
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1 1 l\ 


/i 


1 1 1 
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1 1 1 
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/ill 
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2'2'2'2/\2'2'2'2 


2' 2'2'2/ 
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2'2'2 
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\2'2'2' 
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1 1 1 l\ /l 1 1 1 














2'"2' 2'~2/ \2'~2' 


2' 
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5 

















(158) 



1^14%) = 



1 



tt;t + nn - nn - nn +V2 nn -V2 nn + nn 



+ nn-nn-nn ■ (159) 

Now again for total spin 1, we have 

1^13^24) = ^(|i>i) + |i>o> + |i,-i>), 

where P stands for a pseudoscalar meson. Inserting the completeness expression for an identity / the 
above equation may be written as 

1^3^24)- 4^(10,0, 1,1)(0,0,1,1|)|1,1) + ^(|0,0,1,0)(0,0,1,0|)|1,0) + 



^/3 
1 / 
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(|0,0,1,-1)(0,0,1,~1|)|1,-1). 



Using the Clebsch-Gordan coefficients it becomes 



1^13^24) = ^(|0,0>|1,1> + |0,0)|1,0) + |0,0)|1,-1) 



Using the identity 



1 1 1 l\ /l 1 1 1 




1 1 1 l\ /l 1 1 1 




1 1 1 l\ /l 1 1 1 


2'2'2'2/\2'2'2'2 




2' 2'2'2/\2' 2'2'2 




2'2'2' 2/\2'2'2' 2 



1 _1 1 _1 

2'^2' 2'~2 



1111 
2'^2' 2'"2 



and using the Clebsch-Gordan coefficients we get 



1^4.^2.)-;^^ 



nn - nn +^ ( nn + nn - nn - nn ) 

+ nn - nn 



(160) 



In the same way by inserting the completeness expressions for identity / and using the corresponding 
Clebsch-Gordan coefficients we can have 



1^13^24) = 



nn - nn +^ ( tt;; - nn + nn - nn ) 



1 - 

-^ I •4-"4' ■Y'-Y'-Y' I 



(161) 
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IA.%)^^ 



nn - ;ttt +^ ( tn; + nn - ;nt - ;;tt ) 
+ nil - lilt 



(162) 



1 

4" I YY YY I 4" 



I^i4^23> = ^ [ mt - nn +^ ( m; - nn + nn - nn ) 



V2' 



(163) 



Similarly for total spin we have the fohowing possible spin states for a meson meson molecule 

1 



11^3^24) = 



2 nn +2 nn - nn - nn - nn - nn 



(164) 



1^4%) 



2 nn +2 ;tt; - tt;; - nn - nn - nn 



(165) 



1^13^24) = 2 



nn - nn - nn + nn 



(166) 



I °14^23) ^ 



1 r 



nn - nn - nn + nn 



Overlap of these spin states from eq. (|158p to eq. (J167I) are shown in table [T] and table [51 



For Spin 1 
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(167) 



Table 1: Spin overlaps for spin-1, Meson-meson systems. Here P stands for pseudoscalar and V for vector 
mesons. 



For Spin 
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Table 2: Spin overlaps for spin-0, Meson-meson systems. Here P stands for pseudoscalar and V for vector 
mesons (These agree with eqs. D5 to D8 of appendix D of ref. [25]). 
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Appendix B: The colour basis 



The two meson-like base states {|lj^3l24)c,|l 14 123)0} in terms of orthonormal basis {|3i23B|;)c,|6i26g|;)c} 
are 



|ll3l24>c = -^|3l233|)c+ Yg|6i263j) 



(168) 



|ll4l23)c - -^|3l23 



34/ 



^16126 



34/ 



(169) 



The matrix elements of the identity operator in the new basis {|li3l24)cj|li4l23)c} are shown in 
eq. (PO]) . taking f — I. The matrix elements of colour Casmir operator Fj.Fj in the new colour basis 
{|l;^3l24)c,|li4l23)c} are (see C13 to C15 of appendix- C of ref. [55]) 
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